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The  basic  quantities  usually  determined  during  the  process  of  ignition 
are  the  point  of  ignition  and  the  ignition  delay  period.  We  will  seek  these 
quantities  for  a  drop  of  radius  r0,  placed  in  a  medium  with  temperature  Too 
under  the  following  assumptions:  1)  combustion  of  fuel  vapors  does  not  occur 
until  the  moment  of  ignition  (similar  to  the  assumption  of  0,  M.  Todes  in  the 
unstable  theory  of  thermal  ignition);  2)  Nu=2,  which  is  generally  incorrect, 

since  at  the  initial  moment  in  time  Nu=“>,  after  which  it  rapidly  approaches  its 
stable  value  of  2;  3)  the  times  of  thermal  and  concentration  reaction  are 

infinitesimal  (that  is  the  problem  is  solved  in  the  quasistable  statement); 

4)  there  is  similarity  between  the  fields  of  excess  concentration  and  tempera¬ 
ture  (i.  e.,  the  Lewis-Semenov  criterion  is  equal  to  one);  5)  all  heat  physi¬ 
cal  quantities  are  independent  of  temperature. 

The  quantity  of  heat  used  to  heat  the  drop  is  significantly  greater  than 
the  quantity  of  heat  used  for  evaporation  as  long  as  the  drop  is  far  from  the 
point  of  equilibrium  evaporation.  Therefore,  heating  of  the  drop  can  be  analyzed 
with  some  approximation  without  considering  evaporation.  Assuming  that  the 
temperature  of  the  drop  is  constant  through  its  radius  but  variable  with  time, 
we  produce,  considering  the  assumptions  made, 


Tt 
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Here  T0  is  the  initial  temperature  of  the  drop,  Tj  is  the  (instantaneous)  tem¬ 
perature  of  the  drop,  k  is  the  heat  conductivity  coefiicient,  Cp  and  p  are  the 
heat  capacity  and  density  of  the  drop,  t  is  time.  In  order  to  produce  the  dis¬ 
tribution  of  temperature  in  space,  we  shall  use  the  fundamental  solution  of 
the  stable  heat  conductivity  equation  for  the  spherical  case: 
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Considering  equations  (1)  and  (2),  we  produce,  in  the  quasistable  approximation 


where 
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According  to  the  fourth  assumption  made  above 


r-T,  P»-r± 
1*—  T,  P,»  ‘ 


Here  P,  is  the  partial  pressure  of  fuel  vapors,  a  function  of  the  coordinates 
and  time,  while  Pjq  is  the  partial  pressure  of  the  fuel  vapors  on  the  surface  of 
the  drop,  determined  by  the  saturated  vapor  pressure  at  temperature  Tj.  There¬ 
fore,  the  following  expression  can  be  produced  for  the  distribution  of  fuel  vapor 
pressure  in  space  at  any  moment  in  time: 
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Here  H  is  the  molar  heat  of  evaporation  of  the  fuel,  a  function  of  temperature 
and  pressure,  P  is  the  total  pressure,  T2  is  the  boiling  point  of  the  liquid. 


Ignition  will  occur  at  the  point  of  the  local  maximum  of  the  rate  of  the 
chemical  reaction.  In  all  of  the  following,  we  will  analyze  the  second  order 
chemical  reaction.  Actually,  at  infinity  the  temperature  is  great,  but  the 
fuel  vapor  pressure  is  insufficient.  The  level  of  temperature  near  the  surface 
of  the  drop  is  insufficient.  Furthermore,  the  greatest  number  of  collisions 
leading  to  formation  of  the  end  products  of  the  chemical  reaction  will  occur 
with  the  stoichiometric  ratio.  Therefore,  we  can  write  a  system  of  two  equations 
defining  the  ignition  of  evaporating  liquid  fuel: 
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where  P2  is  the  partial  pressure  of  the  oxidizer  (and  inert  gases) . 

$  is  the  stoichiometric  coefficient, 

the  rate  of  the  chemical  reaction. 

For  a  second  order  reaction,  we  can  approximately  produce,  in  place  of 
what  (1.5)  and  (1.6) 
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where  E  is  the  activation  energy  of  the  chemical  reaction. 


Substituting  the  distributions  of  temperature  and  pressure  found  earlier 
into  (1.7),  we  can  find  the  ratio  of  the  radius  of  ignition  to  the  radius  of 
the  drop  in  the  t-K)  approximation: 
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Following  certain  calculations,  (1.3)  leads  us  to  the  expression 
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under  the  condition  that  _ °  >>  2.  From  this  we  can  find  the  ignition  delay 
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It  follows  from  (1.10)  that  as  temperature  T,*  is  increased,  the  ignition  delay 
period  decreases  and  as  T*,-*®,  t*  0.  As  the  pressure  in  the  system  increases, 
the  magnitude  decreases,  which  in  turn  decreases  t*.  Thus,  the  relationship 
with  pressure  is  concealed  in  the  heat  of  evaporation.  Formula  (1.10)  can  be 
significantly  simplified  if  we  assume  that 


RTi  i 
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Actually,  for  octane,  for  example,  T2—  400°K,  H~10&  cal/mol,  <J>r,l,  &= 

20000  cal/mol.  Then  where  $0  -x  1  and  Tw  =  1100°K  we  have  RT^  ln  F^q  _ 

«  1.  In  this  case  H  RT 
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If  we  assume  T  =300°K,  we  find  using  (1.11)  for  octane  with  2ro=0.1  cm,  c„=0.5 
cal  *g*deg,  p.=0. 8  g/cm3,  k=8.10-5  cal/cm*sec*deg  t*=0.44  sec  and  rj/r0  =0.6  . 

We  can  see  from  the  expressions  produced  that  the  ignition  delay  period  depends 
strongly  on  droplet  size.  The  dependence  on  temperature  is  weaker  than  that 
usually  observed  experimentally,  which  may  be  related  to  the  assumption  of  in¬ 
dependence  of  Cp  and  k  on  temperature. 

We  note  that  in  the  quasistable  approximate  analysis  presented  above,  the 
results  produced  are  independent  of  the  preexponential  factor  of  the  chemical 
reaction.  This  is  explained  by  the  use  of  conditions  (1.5)  and  (1.6). 


For  comparison  with  the  approximate  solution  produced,  the  system  of 
equations  of  diffusion  and  heat  conductivity  was  numerically  calculated  in 


the  form 
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with  the  initial  conditions 


T(r>r„  0)  =  r„  | 
c,(r.  0)  —  0  j 


and  boundary  conditions 


T(os,  t)  —  Tm 
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Here  cj  is  the  fuel  concentration,  c10  is  the  concentration  of  fuel  of  the 
surface  of  a  drop,  Co  is  the  concentration  of  oxidizer,  k0  is  the  preexponential 
factor,  q  is  the  heat  value  of  the  fuel. 


The  following  values  of  dimensional  quantities  were  selected  for  calculation: 
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cp=0.5 
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p  =1.2*10"3  2  -  .  1A3  cal  ;  p  ■  1  atm. 
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Calc..1.,  ju  was  performed  by  the  grids  method  using  the  M-20  computer. 

The  calculations  showed  that  the  increase  in  temperature  begins  at  t=1.385  sec 
at  9  radii  from  the  center  of  the  drop.  Then  the  temperature  begins  to  increase 
significantly,  and  the  temperature  maximum  approaches  the  drop. 

If  we  compare  the  values  of  t*  and  rj  produced  with  the  values  calculated 
according  to  the  approximate  formula-.  '1 .9)  and  H.ll),  we  producer 


£1  *  10. a 

Co 

Tw  s  t  sec. 

The  correspondence  between  approximate  and  precise  calculation  must  be  con¬ 
sidered  quite  satisfactory. 
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